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Gamma- Ray Bursts in 1.8 < z < 5.6 Suggest that the Time 
Variation of the Dark Energy is Small 
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ABSTRACT 

We calibrated the peak energy-peak luminosity relation of GRBs (so called Yonetoku 
relation) using 33 events with the redshift z < 1.62 without assuming any cosmological 
models. The luminosity distances to GRBs arc estimated from those of large amount of 
Type la supernovae with z < 1.755. This calibrated Yonetoku relation can be used as 
a new cosmic distance ladder toward higher redshifts. We determined the luminosity 
distances of 30 GRBs in 1.8 < z < 5.6 using the calibrated relation and plotted the 
likelihood contour in (Q mi Q\) plane. We obtained (fl m ,Q\) — (0.37^']^, 0.63^o'i4) 
for a flat universe. Since our method is free from the circularity problem, we can 
say that our universe in 1.8 < z < 5.6 is compatible with the so called concordance 
cosmological model derived for z < 1.8. This suggests that the time variation of the 
dark energy is small or zero up to z ~ 6. 

Key words: gamma rays: bursts — gamma rays: observation — cosmology: cosmo- 
logical parameters 



1 INTRODUCTION 

There are several distance indicators to determine the cos- 
mological luminosity distance cLl{z). If a certain distance 
indicator is calibrated without any cosmological models, the 
indicator can be used to determine the cosm ological param - 
eters such as fi m , Qa and w = p/p. In 1993. IPhillipsl (|l993l ) 
discovered that the absolute magnitude of Type la super- 
novae at the peak is strongly correlated to the decline rate 
of the lightcurve after the maximum epoch. This correlation 
is rigidly calibrated using other distance indicators such as 
Tully-Fisher and Faber- Jackson relations which are free from 
any cosmological models. If the correlation does not depend 
on the redshift z, we can determine cLl{z) only from the ob- 
served maximum flux and the redshift of the host galaxy. 
Since the luminosity distance depends on the cosmological 
parameters, we can determine them from the estimated lu- 
minosity distances for high redshift Type la supernovae. 

Using the calibrated correlation of Type la supernovae, 
the existence o f the d a rk energy is s t rongl y suggested first by 
ISchmidt etal] i|l99Sf ); iRiess et al.1 (Il998l h iPerlmutter et al.l 
(|l999h . Thanks to the latest large number (41) of observa- 
tions of Type la supernovae, the data with z ^ 1.755 favor 
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(Q m j ^a) = (0.27,0.73) for a flat cosmology (|Riess et al.l 
l2007t ) , which is usually called as the concordance cosmologi- 
cal model. However the most distant Type la supernova ever 
observed is at z — 1.755 so that we need either further Type 
la supernovae or other distance indicators to know the prop- 
erty of the dark energy beyond z > 1.8, while the anisotropy 
of the cosmic microwav e background (CMB ) gives us the in- 
formation at z = 1089 l|Spergel et alj|2007h . 

One of the possible other distance indicators is Gamma- 
ray bursts (GRBs) whose maximum redshift observed is 
higher than that of Type la supernovae. At present the 
most distant GRB is at z = 6.3 confirmed by the spec- 
trosc opic observation with Subaru telescope (|Kawai et al.l 
2006). Since GRBs are known as the most violent and 
brightest explosion in the universe, they might be a 
possible good distance indicator beyond z > 1.8. For 
this purpose, in part, several distance indicators have 
been proposed so far ( Fcnimore & Ramirez- Ruiz .1 
Norris et al. ll2000l:lAmati et al. Il2002l: I Yonetoku et al.l I: 



2000; 



2004: 



Ghirlanda et al.ll2004al : lLiang fc Zhang.ll2005l : iFirmani et al.l 
2006). Using these distance indicators, the cosmological pa- 



rameter s are independently est i mated by seve r al au thors 
such as iGhirlanda et all i|2004bl h IFirmani et~al1 (|2006h and 
ISchaefer et al l (120071 ). 

However the fundamental difficulty exists when we ap- 
ply these distance indicators of GRBs to determine the cos- 
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mological parameters. At first, these indicators are estab- 
lished assuming the concordance cosmological model. Af- 
ter that, cosmological parameters are estimated by various 
methods. Therefore the same parameter set would be in- 
evitably obtained. This logic is falling into a circularity prob- 
lem. Therefore, we should calibrate the distance indicators of 
GRBs without the theoretical concordance cosmology such 
as done for the Type la supernovae. Then we will suggest 
that the cosmic distance ladder is extended, and used for 
the measure ment of the cosmo l ogical constant. 

In 2004, lYonetoku et all l |2004 ) used 11 GRBs with 
known redshifts (0.835 < z < 4.5) at that time assum- 
ing a flat cosmology with S7 m = 0.32, Qa = 0.68 and 
Ho = 72 km s _1 Mpc _1 and derived that the peak lumi- 
nosity (L p ) of GRBs correlates with the peak energy of the 
spectrum E p as L p oc Ep (so called Yonetoku relation). The 
chance probability of this correlation is 5 x 10~ 9 . 

In this Letter, using 63 known redshifts of GRBs, we will 
try to determine the cosmological parameters at higher red- 
shift (1.8 < z < 5.6) only with the observable quantities such 
as the peak flux (f Pt obs) and the E Pl0 bs- Using the luminosity 
distances of many Type la supernovae established without 
any cosmological models, we will calibrate and reconstruct 
the Yonetoku relation for nearby 33 GRBs for z < 1.755. Ap- 
plying the calibrated relation to 30 GRBs in 1.8 < z < 5.6, 
we determine their luminosity distances as a function of z. 
Then we will plot the likelihood contour on the (Q m ,QA) 
plane and obtain Q. m = 0.37^;", !T2a = 0.63±£^ for a flat 
universe. This method (logic) is very similar to the calibra- 
tion for Type la supernovae with the Tully-Fisher and the 
Faber- Jackson relations, and it is free from the circularity 
problem. The essential point is that we use the luminosity 
distances of Type la supernovae like Tully-Fisher relation, 
which are well determined by the observations and are free 
from the cosmological models, to calibrate the Yonetoku re- 
lation. 

In §2 we show how to calibrate the Yonetoku relation 
using the Type la supernovae. In §3 we show the likelihood 
contour in (f2 m ,f2A) plane and argue the cosmological pa- 
rameters. §4 is devoted to discussions. Throughout the pa- 
per, we adopt Ho = 70 km s~ 1 Mpc~ 1 



2 CALIBRATING YONETOKU RELATION 

Generally, the spectrum of the prompt emission of GRBs can 
be explained as an exponentially co nnected broken p ower- 
law model, so called Band function l|Band et al.lll993h . We 
adopted this model to the t ime averaged gamma-r ay spectra 
of known redshift samples (jYonetoku et al.l l2004). Then we 
can determine a peak energy, E p , which co rresponds to the 
energy at the maximum flux in vF v spectra. lYonetoku et al.l 
|2004h found a strong correlation between the E v and the 
1-second peak luminosity (L p ). Here, we estimate the L p 
in 1-10 4 keV energy range in the rest frame of GRB. We 
used the data observed by several independent missions and 
instrumen ts, so we inc l uded appropriate k-correction sug- 
gested by iBloom et alj {2001 ) when we estimate the lumi- 
nosity for each event. We used 33 GRBs for z < 1.755 as the 
calibration of the E p -L p relation. 

We found the empirical form ula of the luminosi ty dis- 
tance of Type la supernovae in iRiess et alj \200H ) with 
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Figure 1. The peak luminosity (L p ) and the peak energy (E p ) 
in the spectrum of 33 GRBs with z < 1.62. The linear correlation 
coefficient is 0.9478 and the chance probability is 6.0 X 10 — 17 . 
The solid line is the best fit curve of L p /10 52 erg s — 1 = 1.31 X 
10- 4 [B P (1 + zVlkeV] 1 - 68 while two dashed lines are the curves 
including the systematic error. See text for details. 

0.359 < z < 1.755 as 

= 14.57 x/' 02 + 7.16 xz 1 - 76 . (1) 

Wcm 

This formula agrees with the real data within a relative error 
of 1 %. Here we note that this formula is not unique and the 
other formula is possible. What is important here is that 
we do not assume any cosmological models at this stage but 
simply assume that the Type la supernovae are the standard 
candle for 0.359 < z < 1.755 irrespective of the cosmological 
model. We apply this formula to 33 GRBs with the redshift 
z < 1.62 in our sample to obtain L p — 4TvdL(z) 2 f p>0 bs while 
E p = (l + z)E Pt0 bs- In figure[T] we show the calibrated E p -L p 
relation for 33 GRBs within z < 1.62. The linear correlation 
coefficient is 0.9478 and the chance probability is 6.12 x 
10 -17 . We tried to find the best fit curve in the form as 

( I5S ^ZT) = (1.31 ± 0.67) xlO"^)— (2) 

In this equation, the error is expressed as the statistical un- 
certainty. However the data distribution has a larger devi- 
ation around the best fit line compared with the expected 
Gaussian distribution. We estimated this systematic devia- 
tion in the normalization as 9.57 x 10 -5 . The solid line is the 
best fit curve and two dashed lines are the curves including 
the systematic error in the normalization. 



3 COSMOLOGICAL PARAMETERS 

At present, there are more than 100 GRBs with known red- 
shift, but the E p has not been measured for all GRBs. More- 
over the calibrated range in figure[T]is 3.3 ^ E p ^ 1333 keV 
in the rest frame of GRBs, so we should use GRBs within 
this range to estimate the luminosity distance of z > 1.755. 
Although we know the redshift of GRB 050904 is z = 6.3, 
its peak energy is reported as E p ~ 3 MeV. This is out 
of our calibration, so we exclude it from our sample. Then, 
we apply the calibrated Yonetoku relation to 30 GRBs in 
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Figure 2. The luminosity distance as a function of the redshift 
measured by the calibrated E p -L p relation. The blue and the red 
points are the luminosity distance of z sC 1.755 and z > 1.755, 
respectively. The data of Type la supernovae are also plotted as 
the black cross points. The uncertainty bar of each point includes 
the systematic dispersion around the best fit line of figure [T] For 
the purpose of easy comparison, a solid line of the concordance 
model is also drawn. 



1.8 < z < 5.6 to determine the luminosity distance as a 
function of z. In figure 2 we show the luminosity distance 
of 33+30 GRBs as a function of the redshift. The blue and 
the red points are the luminosity distance of z ^ 1.755 and 
z > 1.755, respectively. The uncertainty bar of each red 
point includes the systematic dispersion around the best fit 
line of figure [1] 

For each GRB with z = z 1 we have the observed peak 
flux(/p o6s ) in the unit of erg cm _2 s -1 and the observed 
E Pt obs m the unit of 1 keV. Then using the equation [2] the 
luminosity distance can be derived as 



<1l(z 1 ) — 10 24 cm 



1.31 



-[*4 oi , s (i + * i )] 1 - 88/2 - 



(3) 



In the ACDM-cosmology with fi*. 
nosity distance is given by 



A 








Figure 3. The contour of the likelihood A\ 2 (see equation[6]l for 
the luminosity distances of 30 GRBs in 1.8 < z < 5.6. The sig- 
nificance levels of 68%, 95% and 99% are also shown on the same 
panel. Compared with the case of Type la supernovae, the shape 
of the contour stands vertical since the mean redshift is much 
larger. The most likelihood value of cosmological parameters are 
(n m ,n A ) = (0.251q 2 4, 1.25l5'25) while for a flat cosmology prior 
they are (« m , « A ) = (0.37!^, 0.63±2;"). 



d\ (z, f2 m , Q\) 



where F(z) = 



Qa — 1, the lumi- 
sin(V£h.F(z)) if Q k > 



(0.25±o; 2 4, 1.25+° l 5 ) while for a flat cosmology prior that 
is (fi m ,fi A ) = (0.37tg:n,0.63±|^ 4 ) with 1 a uncertainty. 

In figure [4j we show the same contour as figure [3] but 
for the luminosity distances of 16 GRBs in 3 < z < 5.6. 
The shape of the contour stands more vertical than fig. [3] 
which should be so. Without any prior the most likeli- 
hood values of cosmological parameters are (fi m ,fi A ) = 
(0.33to jg, T14t?; 2 4) while for a flat cosmology prior they 
are (« m ,Q A ) = (6.49± ;||, 0.53^; 22 ). We see for a fiat cos- 
mology the value of (fi m , S7 A ) for higher redshift samples is 
compatible with that for whole samples if we take the error 
into account. 



: smh(y^nTF( 2 )) if Q k < 
if fl k = 

dz 
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4 DISCUSSIONS 



/() v/n m (l+2) 3 +fiA-«fc(l + z) 2 

We define a likelihood function as 



30 

» 2 _^ Aogd L (z l )-i og d t L h (z\n m ,n A ) 2 

X ^ ( Ad L {z*) ' 

i — 1 



2 

Xbest • 



(6) 



Here Xbest means the chi-square value for the best fit pa- 
rameter set of Q m and Qa- In figure [3] we show the con- 
tour of the likelihood A\ 2 for the luminosity distances 
of 30 GRBs in 1.8 < z < 5.6. Compared with the case 
of Type la supernovae, the shape of the probability con- 
tour stands vertical since the luminosity distance strongly 
depends on fi m rather than f2 A for higher redshift sam- 
ples. This is clear from the functional form of F(z). With- 
out any prior the most likelihood value is (f2 m ,fl A ) = 



(5) In this Letter we extended the cosmic distance ladder using 
GRBs calibrated by Type la supernovae, and argued the cos- 
mological parameters using only the GRBs in 1.8 < z < 5.6 
except for a flat cosmology prior. Since our region of 1.8 < 
2 < 5.6 has not been explored yet, this is the first report 
to estimate the cosmological parameters up to z — 5.6. It is 
possible to combine our data with (1) Type la supernova, 
(2) CMB, (3) Baryon Acoustic Oscillation (4) the large scale 
structure measurement and (5)weak gravitational lensing, to 
constrain the w(z) = p/p parameter of the equation of state 
for the dark energy (Tsutsui et al.). 

The calibrated Yonetoku relation has a large disper- 
sion in the normalization which is mainly caused by the 
systematic error. Therefore, currently, the measurement of 
the luminosity distance is not so accurate compared with 
the other distance indicators such as Type la supernova. It 



4 Y. Kodama et al. 




Figure 4. The same as fig. [3] but for the luminosity distances of 
16 GRBs in 3 < z < 5.6. The shape of the contour stands more 
vertical than fig. [2] The most likelihood values are (Q. m ,&<A.) = 
(0.33^ '2g , 1.14^j while for a flat cosmology prior they are 
(Q m ,n A ) = (0.49t°;g,0.53l°;^). 



may be difficult to discuss the detailed time history of the 
cosmological parameters yet. However, if this deviation is 
the intrinsic property of GRBs, we will be able to discover 
the hidden physical quantities, or distinguish a possible sub- 
class from entire population of GRBs like Type la supernova. 
These improvements in the cosmic distance ladder will lead 
us to explore the deep space with better accuracy in near 
future. 

As for the Amati-relation (|Amati et al . 2002) which 
relates the total isotropic energy to E v , possible se lec- 
tion bias effects and e volutio n effects are cl a imed l|Li.l 
120071 ; iButler et all I2007T ) while IWillingale et~ai] (|2007T ) ar- 
gues against such effects. It is important to check possible 
selec tion bias and evoluti on effects in the Yonetoku relation 
also. iTanabe et al.l (|2008l ) examined the evolution effect as 
well as the observational selection bias assuming the concor- 
dance cosmology, and found that they are quite small. We 
think that the selection bias and evolution effects are not 
larger than the systematic uncertainty in the normalization 
of the calibrated Yonetoku relation, although we agree that 
further examinations are required. 

Theoretical models of dar k energy are reviewed, for 
examp le, in a recent paper by iFrieman. Turner fc Hutereij 
(2008). In some models such as scalar field models, dark 
energy looks like the cosmological constant for low redshift 
2 < 2 but not for high 2 > 2. Therefore the estimate of 
the cosmological parameters for 2 > 2 is important either 
to refute or confirm such models. Our present result for 
1.8 < 2 < 5.6 gives (On,Ov) = (0.37±g;^, 0.63±°;^) for a 
flat cosmology prior while Type la supernova for 2 < 1.755 



does (O ro ,fi A ) = (0.27,0.73) (|Riess et al l (|2007n ), which 
means that (Sl m ,OA) for 2 < 1.8 and for 1.8 < 2 < 5.6 
are the same within the 1 sigma statistical error. This sug- 
gests that (fi m , Qa) did not change so much from 2 = 5.6 to 
2 = so that such a scalar fieldmodel should look like the 
cosmological constant up to 2 ~ 6 at least. 
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